rr library(dplyr)


Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union

rr wolfrum <- readRDS(‘/projects/timshel/sc-scheele_lab_adipose_fluidigm_c1/data-wolfrum/wolfrum.compute.seurat_obj.rds’)

What do the clusters represent?

rr pos_markers_top20 <- markers %>% group_by(cluster) %>% top_n(n=6, wt=avg_logFC) %>% unite(col=‘cluster_gene’, c(cluster, gene), remove=F)

neg_markers_top20 <- markers %>% group_by(cluster) %>% top_n(n=6, wt=desc(avg_logFC)) %>% unite(col=‘cluster_gene’, c(cluster, gene), remove=F)

rr plots <- FeaturePlot(wolfrum, features=as.vector(pos_markers_top20\(gene[1:6]), combine = FALSE, pt.size=1) for(i in 1:length(plots)) { plots[[i]] <- plots[[i]] + ggtitle(pos_markers_top20\)cluster_gene[i]) + NoAxes() }

rr plot_grid(plotlist=plots, ncol=3)

U and L branch markers.

rr FeaturePlot(wolfrum, features=c(‘UCP2’, ‘FABP5’, ‘G0S2’, ‘FABP4’, ‘ADIPOQ’, ‘ADIRF’, ‘CD36’, ‘PLIN4’, ‘APOD’, ‘MGP’, ‘DCN’, ‘CTGF’, ‘IGF2’, ‘CCDC80’, ‘PLAC9’, ‘THBS1’), cols=c(‘gray’, ‘blue’), ncol = 2, reduction=‘umap’, pt.size=1)

Preadipocytes

rr FeaturePlot(wolfrum, features=c(‘TMSB4X’), cols=c(‘gray’, ‘blue’), ncol = 2, reduction=‘umap’, pt.size=1)

Monocle

Ran Monocle on the Wolfrum subset with the same gene list that was used to build the 10x-180831 trajectory.

monocle_results <- readRDS('../output/monocle/wolfrum/wolfrum_subset.compute.seurat_obj.rds-monocle-monocle_genelist_T1T2T3_T4T5_res.1.5')

  return(df)
Error: no function to return from, jumping to top level
states.prediction
     n Metabolic  ECM Progenitor
1 3242      0.96 0.01       0.03
2 2385      0.09 0.84       0.07
3 1605      0.64 0.09       0.28
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKCmBgYHtyfQpsaWJyYXJ5KFNldXJhdCkKbGlicmFyeShjb3dwbG90KQpsaWJyYXJ5KGRwbHlyKQpsaWJyYXJ5KHRpZHlyKQpgYGAKCmBgYHtyfQp3b2xmcnVtIDwtIHJlYWRSRFMoJy9wcm9qZWN0cy90aW1zaGVsL3NjLXNjaGVlbGVfbGFiX2FkaXBvc2VfZmx1aWRpZ21fYzEvZGF0YS13b2xmcnVtL3dvbGZydW0uY29tcHV0ZS5zZXVyYXRfb2JqLnJkcycpCmBgYAoKYGBge3IgZmlnLmhlaWdodCA9IDYsIGZpZy53aWR0aCA9IDEyLCBmaWcuYWxpZ24gPSAiY2VudGVyIn0KcGxvdF9ncmlkKAogIFVNQVBQbG90KHdvbGZydW0sIGdyb3VwLmJ5PSdvcmlnLmlkZW50JywgbGFiZWw9VCksCiAgVU1BUFBsb3Qod29sZnJ1bSwgZ3JvdXAuYnk9J3NldXJhdF9jbHVzdGVycycsIGxhYmVsPVQpCikKYGBgCgpXaGF0IGRvIHRoZSBjbHVzdGVycyByZXByZXNlbnQ/CgpgYGB7cn0KbWFya2VycyA8LSByZWFkLnRhYmxlKCcuLi9vdXRwdXQvbWFya2VyZ2VuZXMvd29sZnJ1bS9tYXJrZXJzX3dvbGZydW0uY29tcHV0ZS5zZXVyYXRfb2JqLnJkc19zZXVyYXRfY2x1c3RlcnNfbmVnYmlub20nLCBzZXA9J1x0JywgaGVhZGVyPVQpCmBgYAoKYGBge3J9CnBvc19tYXJrZXJzX3RvcDIwIDwtIG1hcmtlcnMgJT4lIAogIGdyb3VwX2J5KGNsdXN0ZXIpICU+JSAKICB0b3BfbihuPTYsIHd0PWF2Z19sb2dGQykgJT4lIAogIHVuaXRlKGNvbD0nY2x1c3Rlcl9nZW5lJywgYyhjbHVzdGVyLCBnZW5lKSwgcmVtb3ZlPUYpCiAgCm5lZ19tYXJrZXJzX3RvcDIwIDwtIG1hcmtlcnMgJT4lIAogIGdyb3VwX2J5KGNsdXN0ZXIpICU+JSAKICB0b3BfbihuPTYsIHd0PWRlc2MoYXZnX2xvZ0ZDKSkgJT4lCiAgdW5pdGUoY29sPSdjbHVzdGVyX2dlbmUnLCBjKGNsdXN0ZXIsIGdlbmUpLCByZW1vdmU9RikKYGBgCgpgYGB7ciwgZmlnLnNob3c9J2hpZGUnLCByZXN1bHRzPSdoaWRlJywgbWVzc2FnZT1GfQpwbG90cyA8LSBGZWF0dXJlUGxvdCh3b2xmcnVtLCBmZWF0dXJlcz1hcy52ZWN0b3IocG9zX21hcmtlcnNfdG9wMjAkZ2VuZVsxOjZdKSwgY29tYmluZSA9IEZBTFNFLCBwdC5zaXplPTEpCgpmb3IoaSBpbiAxOmxlbmd0aChwbG90cykpIHsKICBwbG90c1tbaV1dIDwtIHBsb3RzW1tpXV0gKyBnZ3RpdGxlKHBvc19tYXJrZXJzX3RvcDIwJGNsdXN0ZXJfZ2VuZVtpXSkgKyBOb0F4ZXMoKQp9CmBgYAoKYGBge3IsIGZpZy5oZWlnaHQgPSAxMTIsIGZpZy53aWR0aCA9IDEwLCBmaWcuYWxpZ24gPSAiY2VudGVyIn0KI2hlaWdodCAxIHJvdyA9IDMuIAojMjcgY2x1c3RlcnMuIDI3KjIgPSA1NCByb3dzLiBoZWlnaHQgPSA1NCozID0gMTEyIApwbG90X2dyaWQocGxvdGxpc3Q9cGxvdHMsIG5jb2w9MykKYGBgCgpVIGFuZCBMIGJyYW5jaCBtYXJrZXJzLiAKCmBgYHtyIGZpZy5oZWlnaHQgPSAzMCwgZmlnLndpZHRoID0gOSwgZmlnLmFsaWduID0gImNlbnRlciJ9CkZlYXR1cmVQbG90KHdvbGZydW0sIGZlYXR1cmVzPWMoJ1VDUDInLCAnRkFCUDUnLCAnRzBTMicsICdGQUJQNCcsICdBRElQT1EnLCAnQURJUkYnLCAnQ0QzNicsICdQTElONCcsICdBUE9EJywgJ01HUCcsICdEQ04nLCAnQ1RHRicsICdJR0YyJywgJ0NDREM4MCcsICdQTEFDOScsICdUSEJTMScpLCBjb2xzPWMoJ2dyYXknLCAnYmx1ZScpLCBuY29sID0gMiwgcmVkdWN0aW9uPSd1bWFwJywgcHQuc2l6ZT0xKQpgYGAKClByZWFkaXBvY3l0ZXMKCmBgYHtyIGZpZy5oZWlnaHQgPSAzMCwgZmlnLndpZHRoID0gMTIsIGZpZy5hbGlnbiA9ICJjZW50ZXIifQpGZWF0dXJlUGxvdCh3b2xmcnVtLCBmZWF0dXJlcz1jKCdUTVNCNFgnKSwgY29scz1jKCdncmF5JywgJ2JsdWUnKSwgbmNvbCA9IDIsIHJlZHVjdGlvbj0ndW1hcCcsIHB0LnNpemU9MSkKYGBgCgojTW9ub2NsZQoKUmFuIE1vbm9jbGUgb24gdGhlIFdvbGZydW0gc3Vic2V0IHdpdGggdGhlIHNhbWUgZ2VuZSBsaXN0IHRoYXQgd2FzIHVzZWQgdG8gYnVpbGQgdGhlIDEweC0xODA4MzEgdHJhamVjdG9yeS4gCgpgYGB7cn0KbW9ub2NsZV9yZXN1bHRzIDwtIHJlYWRSRFMoJy4uL291dHB1dC9tb25vY2xlL3dvbGZydW0vd29sZnJ1bV9zdWJzZXQuY29tcHV0ZS5zZXVyYXRfb2JqLnJkcy1tb25vY2xlLW1vbm9jbGVfZ2VuZWxpc3RfVDFUMlQzX1Q0VDVfcmVzLjEuNScpCgpgYGAKCmBgYHtyfQpwbG90X2dyaWQoCiAgcGxvdF9jZWxsX3RyYWplY3RvcnkobW9ub2NsZV9yZXN1bHRzLCBjb2xvcl9ieT0nU3RhdGUnKSwKICBwbG90X2NlbGxfdHJhamVjdG9yeShtb25vY2xlX3Jlc3VsdHMsIGNvbG9yX2J5PSdwcmVkLmJyYW5jaC5pZCcpCikKYGBgCgpgYGB7cn0Kd29sZnJ1bV9zdWJzZXQgPC0gQWRkTWV0YURhdGEod29sZnJ1bV9zdWJzZXQsIHBEYXRhKG1vbm9jbGVfcmVzdWx0cylbJ1N0YXRlJ10pCmBgYAoKCmBgYHtyfQpnZXRfcmF0aW9zIDwtIGZ1bmN0aW9uKGRhdGEsIGNvbDEsIGNvbDIpewogIHN0YXRlcyA8LSB1bmlxdWUoZGF0YUBtZXRhLmRhdGFbLGNvbDFdKQogIHZhbHVlcyA8LSB1bmlxdWUoZGF0YUBtZXRhLmRhdGFbLGNvbDJdKQogIGRmIDwtIGFzLmRhdGEuZnJhbWUobWF0cml4KG5jb2w9bGVuZ3RoKHZhbHVlcykrMSwgbnJvdz1sZW5ndGgoc3RhdGVzKSkpCiAgY29sbmFtZXMoZGYpIDwtIGMoJ24nLCB2YWx1ZXMpCiAgcm93bmFtZXMoZGYpIDwtIHN0YXRlcwogIGZvciAoc3RhdGUgaW4gc3RhdGVzKXsKICAgIG5fc3RhdGUgPSBsZW5ndGgod2hpY2goZGF0YUBtZXRhLmRhdGFbY29sMV0gPT0gc3RhdGUpKQogICAgZGZbc3RhdGUsICduJ10gPC0gbl9zdGF0ZQogICAgZm9yICh2YWx1ZSBpbiB2YWx1ZXMpewogICAgCW5fc3RhdGVfdmFsdWUgPC0gbGVuZ3RoKHdoaWNoKGRhdGFAbWV0YS5kYXRhW2NvbDFdID09IHN0YXRlICYgZGF0YUBtZXRhLmRhdGFbY29sMl0gPT0gdmFsdWUpKQoJICAgIHBlcmNfc3RhdGVfdmFsdWUgPC0gbl9zdGF0ZV92YWx1ZSAvIG5fc3RhdGUKCSAgICBkZltzdGF0ZSwgdmFsdWVdIDwtIHJvdW5kKHBlcmNfc3RhdGVfdmFsdWUsIDIpCiAgICB9CiAgfQogIHJldHVybihkZikKfQpgYGAKCmBgYHtyfQpzdGF0ZXMucHJlZGljdGlvbiA8LSBnZXRfcmF0aW9zKGRhdGE9d29sZnJ1bV9zdWJzZXQsIGNvbDE9J1N0YXRlJywgY29sMj0ncHJlZC5icmFuY2guaWQnKQpzdGF0ZXMucHJlZGljdGlvbgpgYGAKCmBgYHtyfQoKYGBgCgoK